function [rho, Sigma_v] = drawSerialCorrelations_new_Simple(rho,idiosync,Sigma_v,rho_priors)
%UNTITLED6 Summary of this function goes here
%   Detailed explanation goes here
            
            n = size(idiosync,2);
            
            
            idiosync_stand = idiosync./(sqrt(Sigma_v)');  
            sigma2 = 1; 
            
            
            
            for i=1:n                                  
              
                u = idiosync_stand(:,i);
                ulags = mlag2(u,2);
                
                % Adjust for Missing Observations
                
                keep=isfinite(sum([u ulags],2));
                ulags = ulags(keep,:);
                u = u(keep);
                
                % Draw rho but  ensure stationarity
                
                
                
                
                chck=-1;    
                while chck<0               
                [RHO,sigma2] = bayes_reg(u,ulags,rho(i,:),sigma2,rho_priors,0);        
                temp=[RHO';eye(2-1,2)];
                ee=max(abs(eig(temp)));
                if ee<=1; chck=1; 
                else
                    disp('Warning: Explosive Root in Idiosyncratic Component')
                    disp(i);
                    accept = rand(1) <= 0.1;
                    if accept == 1; chck=1; end
                    
                end;
                end
                
                rho(i,:) = RHO;
                
            end
            


end

